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plane is found to be more extended than that of mechanical instability. However, 
locations of these boundaries depend sensitively on both the bulk compressibility 
of nuclear matter and the density dependence of the symmetry energy. 
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I. INTRODUCTION 


During the last decade, there have been extensive studies on multifragmentation in nu¬ 
clear reactions, which is characterized by the emission of several intermediate mass fragments 
with a power law distribution. It has been suggested that this phenomenon is due to various 
mechanical instabilities of the nuclear matter, such as the volume instability, the surface 
instability of the Rayleigh kind, and the Coulomb instability. A recent review on nuclear 
multifragmentation can be found in ref. |IJ. With radioactive ion beams, which have been 
successfully developed recently, nuclear matter with a high isospin asymmetry can be created 
transiently during reactions induced by neutron-rich nuclei (e.g., @.0). This thus offers the 
possibility to study also the chemical instability associated with the isospin asymmetry of the 
nuclear matter. Indeed, recent experiments on nuclear reactions involving different ratios of 
the neutron to proton numbers have shown that products from nuclear multifragmentation 
depends strongly on the isospin asymmetry of the colliding nuclei 0-^|] . 

It is well-known that nuclear matter is not thermodynamically stable at all densities 
p, temperatures T, and neutron excess 5 = ( p n — p p )/p■ The necessary and sufficient 
conditions for the existence of a stable isospin-asymmetric nuclear matter can be expressed 
by the following inequalities (e.g., [i7j-0): 



( 1 ) 


( 2 ) 


( 3 ) 


where E, P, and p, n are the energy per nucleon, pressure and neutron chemical potential, 
respectively. The first condition is required by the thermodynamic stability and is always 
satisfied by any reasonable nuclear equation of state; the second condition ensures the me¬ 
chanical stability against the growth of density fluctuations; and the last one protects the 
chemical or diffusive stability (DS) against the development of neutronization. The last two 
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conditions can be violated in some regions of the (p, T, 5) configuration space, leading to 
mechanical and/or chemical instabilities. 

The boundaries of mechanical and chemical instabilities in the (p, T, 5) space depend on 
the equation of state (EOS) of asymmetric matter. In the Relativistic Mean-Field (RMF) 
theory, which has a linear density dependence in the symmetry energy, recent studies indicate 
that the diffusive spinodal associated with the chemical instability of asymmetric nuclear 
matter encloses more of the configuration space than the isothermal spinodal associated with 
the mechanical instability [|9]. For values of isospin asymmetry, S < 0.4, which are reachable 
in reactions induced by neutron-rich nuclei, the separation between the isothermal and 
diffusive spinodals in P — p plane is significant, and the diffusive instability becomes a more 
relevant spinodal in discussing the properties of asymmetric nuclear matter. 

However, the density dependence of nuclear symmetry energy is poorly known, and 
theoretical predictions vary widely mm- Although it is linear in density in the RMF 
theory, a p 1 / 3 dependence was obtained by Siemens using the Bethe-Goldstone theory for 
asymmetric nuclear matter [|TUJ . More sophisticated calculations by Wiringa et al. using 
the variational many-body theory |T3| give different density dependence depending on the 
nuclear forces used in the calculations. It is thus both interesting and necessary to examine 
how the structure of mechanical and chemical instabilities in asymmetry nuclear matter 
might change when other nuclear equations of state are used. 

The nuclear symmetry energy at densities other than the normal nuclear matter density 
is of interest in both astrophysics and nuclear physics. For example, the cooling rate, chemi¬ 


cal composition and magnetic moment of a neutron star depend critically on this [15,17,18 


For nuclear reactions induced by neutron-rich nuclei, it was recently found that the nuclear 
symmetry energy affects significantly the reaction dynamics, especially the ratio of preequi¬ 
librium neutrons to protons emitted in the early stage of heavy-ion collisions at intermediate 
energies Ell¬ 
in the present paper, we study the chemical and mechanical instabilities of asymmetric 
nuclear matter using various phenomenological equations of state. These EOS’s have been 
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used previously by us in studying the ratio of preequilibrium neutrons to protons in nuclear 
reactions induced by neutron-rich nuclei. In particular, we examine the dependence of the 
boundaries of these instabilities on both the compressibility and the density dependence 
of the symmetry energy. We find that for all equations of state considered here a similar 
conclusion as that based on the RMF is obtained, i.e., the isothermal spinodal is enclosed 
within the diffusive spinodal, and the location of these two spinodals depends sensitively on 
both the bulk compressibility of nuclear matter and the density dependence of the symmetry 
energy. 


II. THERMODYNAMICAL ANALYSIS OF CHEMICAL AND MECHANICAL 

INSTABILITIES 


For asymmetric nuclear matter in thermal equilibrium, the density p q of neutrons (q = n) 
or protons (q = p ) is 

i r°° 

Pq = — / k 2 f q (k)dk, (4) 

7 + JO 

where 


fg(k) = [exp(e g - p q )/T + 1] 1 


( 5 ) 


is the Fermi distribution function, and e q is the single particle kinetic energy. For modestly 
high temperatures (T > 4MeV), this equation can be inverted analytically to obtain the 
chemical potential |20| [22| 


Pq 


V q + T 


ln(%) + E 

Z 71=1 


71+1 

Un 

n 



( 6 ) 


where 


A t — 


'2tt h 2 
mT 


( 7 ) 


is the thermal wavelength of the nucleon, and b' n s are the inversion coefficients given in refs. 
In eq. (||), V q is the single particle potential energy and can be parameterized as 
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V q (p, S) = a(p/p NM ) + Kp/pnmY + V*(p, S). 


( 8 ) 


The parameters a, b and cr are determined by the saturation properties and the compress¬ 
ibility K of symmetric nuclear matter, i.e., 

K + 44.73 


a =-29.81-46.90 A ._ m32 
6 = 23.45 ^ 25578 (MeV), 


(MeV), 


a = 


K - 166.32 
K + 44.73 


(9) 

( 10 ) 

( 11 ) 


211.05 

The isospin-independent term should also contain a momentum-dependent part which 
is important for some dynamical observables, such as the collective flow (e.g., p3|-|25|), but 
is not essential for other observables, such as the ratio of preequilibrium neutrons to pro¬ 
tons []19|. Effects of momentum-dependent interactions on thermal properties were studied 
recently in refs. |2f]j27|] and were found to have only small effects. We thus neglect the 
momentum-dependent interaction and concentrate on investigating effects of the isospin- 
asymmetric potential V a | y . This potential is given by 

Ksy(P> <*) = dw a (p, 8)/dp q , (12) 

where w a (p,5) is the contribution of nuclear interactions to the symmetry energy density, 


i.e., 


and 


w a (p,S) = e a ■ pF(u)5 2 , 


e a = S 0 - (2 2 / 3 - 1 )-E° f . 

5 


(13) 


(14) 


In the above, u = p/po with p 0 !h e normal nuclear matter density; S 0 is the symmetry 
energy S(p) at po and is known to be in the range of 27-36 MeV 0; and Ep is the Fermi 
energy at po- 

For the function F(u ), we consider the following three forms: F(u) = n 2 , u and n 1//2 , 
which are taken form the results of typical microscopic many-body calculations. Then the 
corresponding symmetry potentials are, respectively, 
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(15) 

(16) 


V'"<" ) = ±2e„« 2 a + e a u 2 S\ 
= ±2e 0 ui, 


and 


TAn(p) 

asy 


±2e u it 1/2 5 - ^e a ii 1/2 3 2 , 


(17) 


where “+” and ” are for neutrons and protons, respectively. 

From the chemical potentials determined by eq. (fjj), we obtain the total pressure of the 
system from the Gibbs-Duhem relation, 


dP 

dp 


P 

2 


(1 + 5 ) 


d pn 

dp 


+ (1-5) 


dp p 

dp 


The result can be separated into three parts, i.e., 


(18) 


P — +kin + +0 + Pasy j 


(19) 


where P^ n is the kinetic contribution 


fldn = Tp{ 1 + i y 5„(+^r [(1 + fl" 41 + (1 - <5) 1+ "]}, 

Z n =1 4 

and P 0 is the contribution from the isospin-independent nuclear interaction 


( 20 ) 


p° = +p»(-) 2 +T+-po(-r +1 - 

2 Po 1 + CT Po 


( 21 ) 


P asy is the contribution from isospin-dependent part of the nuclear interaction, and it can 
be written as 


Pasy = 2e a p 0 ( —) 3 <f 2 , 
Po 

Pasy = e a p 0 ( —) 2 <f 2 , 

Po 

and 

Pasy = ^e a p 0 (-j-) 3/2 <5 2 
^ Po 

for F(w) = u 2 , u and n 1//2 , respectively. 


( 22 ) 

(23) 


( 24 ) 
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The condition for mechanical stability follows from eq. fll8P and can thus be determined 
straightforwardly. To evaluate the condition for chemical stability, we use the following 
relations: 


d/-in \ (dp n \ ( dP \ 

d6) Tp -{d6) Tp \dp) T /{d6) T ; 

_ fdpA t —} 1 (—\ 

\ d6 ) T , P V d P / T,5 \ d p) T ,S \ d5 JT,p' 


(25) 


Knowing and (f“) r<5 as functions of p, T and 5, we can then identify regions of 

the configuration space where chemical and/or mechanical instability occur. We first show in 
Fig. 0 and Fig. [| these two quantities as functions of p at various T and 5. In the calculation 
K = 200 MeV and F(u) = u are used. For comparison we also show in the left panel of 
Fig. | results for symmetric nuclear matter. In this case there is no chemical instability and 
the quantity (^f) increases with both density and temperature. To understand this, we 
note that 


^ = ]p™, + P f>+1 )M^r «i+i) n - (i - <*)”). 


' 3 p n 

35 


T,P 


'dV£y' 

35 


+ T 


T,p 


1 +E(«+W^'f (i+^r 1 


1 + 6 


n =1 


and 


(26) 

(27) 


Then, it can be shown that 


linp^o 



0. 


(28) 


lim 5 _ 


3pn 

~35 


= hm 5 _ 


T,P 


3 Pn 
~35 


= 2 e a F(u) + T 


1 + (n + 1)6,, 


T,p 

( -6 tP 


n =1 


(29) 

(30) 


which is always positive as expected for symmetric nuclear matter. Furthermore, it increases 
with both density and temperature. 

Mechanical instability is known to happen at intermediate densities corresponding to 
the mixed phase between the gas and liquid phases at subcritical temperatures. This is 
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clearly demonstrated by plotting as a function of p. It is interesting to examine 

this quantity at the low density limit. From the Gibbs-Duhem relation of eq. (|18|) and the 
expression for p q of eq. (|6|), we obtain 


lim 


'dP s 


p-> o 


= lim 


p^o(T + a—) — T, 
Po 


(31) 


\ 9 P / T,S 

as in a dilute gas. 

We note from eq. ( |25| ) that the quantity is singular along the boundary of 

mechanical instability where (j^) Tg = 0- This singularity is indicated by the spikes in Fig. 
jl] and Fig. Q. However, the height of these spikes between the boundaries of the gas phase 
and the mixed phase is relatively low. This is because at this boundary the last derivative 
in eq. (|25|) is also close to zero, i.e., 


(dP\ 2 

linip^o I ) = - • lim^^nT’^v = 0. 


35 ) 


p—>0-* asy 


T,p 


Although a nuclear system is always chemically stable in the low density limit as 


lim 


p —>o 


' 3\i n 
35 


T 


P,T 


1 + 5 


> 0 , 


(32) 


(33) 


it starts to develop chemical instability at finite density once it becomes neutron-rich. In 
particular, it is seen that the boundary of mechanical instability shrinks while that of chemi¬ 
cal instability expands as the isospin asymmetry 5 increases. For example, at a temperature 
of T=8 MeV the mechanical instability gradually disappears but the chemical instability 
becomes more pronounced as 5 increases from 0 to 0.6. Furthermore, in the mechanically 
unstable region the nuclear system can be chemically stable, but at higher densities chem¬ 
ical instability happens even in the mechanically stable region. This is more clearly shown 
in Fig. [3] where the pressure along the chemical (DS) and mechanical (ITS) spinodals are 
plotted as functions of density. Both the isothermal pressure and separation between the 
two spinodals increase with increasing 5. Therefore, not only the chemical but also the 
mechanical instability is strongly isospin-dependent. 

The isospin dependence of mechanical and chemical instabilities at a fixed temperature is 
studied in more detail in Fig. (| and Fig. || It is seen that the system is both mechanically 
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and chemically stable at low and high densities. At intermediate densities (e.g., p/po = 
0.4, 0.5) the system is mechanically unstable at low 5 (e.g., S < 0.2 and 5 < 0.4 for p/p 0 = 0.4 
and 0.5, respectively) and chemically unstable at intermediate 5 (e.g., 0.2 < 5 < 0.4 for 
p/po = 0.4 and 0.4 < 5 < 0.6 for p/p 0 = 0.5, respectively). The corresponding boundaries 
in the pressure-density plane are shown in Fig. ^|. Again, the diffusive spinodal line is 
more extended than the isothermal spinodal. The isospin dependence of the pressure is also 
shown in this figure. We note that its strong dependence on 5 is due to the significant Pasy 
contribution to the total pressure. 

Both chemical and mechanical instabilities also depend on the isospin-dependent and 
-independent parts of the nuclear equation of state. To demonstrate this, we show in Fig. 
j6] and Fig. i7] the pressure as a function of density p (left panels) and isospin asymmetry S 
(right panels) along the diffusive spinodals (upper windows) and isothermal spinodals (lower 
windows) at a constant temperature of T = 10 MeV by using the three forms of F(u ) and 
a compressibility of 380 and 200 MeV, respectively. From these two figures we observe the 
following interesting features. First, the diffusive spinodals are always more extended than 
the isothermal spinodals. Second, both DS and ITS depend on the form of F(u), and this 
dependence is stronger for the stiff equation of state of K = 380 MeV. Third, from the right 
panels of these two figures one sees that it is both chemically and mechanically favorable 
for the system to be less asymmetric (smaller 6) in the liquid phase than in the gas phase. 
This result is consistent with that based on the energy consideration. Since the equation 
of state for asymmetric nuclear matter contains a S(p)S 2 term, it is therefore energetically 
favorable for it to separate into a liquid phase that is less asymmetric and a gas phase that 
is more asymmetric, rather than into two phases with equal isospin asymmetry. Finally, 
it is seen that there is a maximum isospin asymmetry for both DS and ITS. Along both 
DS and ITS the isospin asymmetry increases to a maximum value and then decreases as 
the density is reduced from the liquid phase towards the gas phase. As shown in both 
figures, the maximum isospin asymmetry is also sensitive to both the isospin-dependent and 
-independent parts of the nuclear equation of state. 
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III. SUMMARY AND DISCUSSIONS 


Using various phenomenological equations of state, we have studied both the chemical 
and mechanical instabilities in hot isospin-asymmetric nuclear matter. Both instabilities 
are found to depend strongly on the isospin asymmetry of the nuclear matter. Also, in 
a chemically and mechanically stable asymmetric nuclear matter a lower neutron excess 
is favored in the liquid phase than in the gas phase. Furthermore, there is a maximum 
neutron excess in the mixed phase at intermediate densities. For all nuclear equations of 
state considered here, the boundary of chemical instability is found to be more extended 
than that of mechanical instability. However, the location of these boundaries in the (T, 
p, <5) space depends sensitively on both the bulk compressibility of nuclear matter and the 
density dependence of the symmetry energy. These results are relevant in understanding the 
significant differences observed in the fragment distributions from collisions involving isospin- 
symmetric and isospin-asymmetric nuclei at intermediate energies [jj-|(|. One scenario for 
nuclear multifragmentation is that the hot system formed in a reaction expands almost 
adiabatically into the mechanical instability region and then disassemble into clusters and 
nucleons due to the growth of density fluctuations (e.g. [|^]). Great efforts have thus been 
devoted during the last decade (e.g., |30|-p4|) to incorporate the effects of density fluctuations 
into dynamical transport models. To describe nuclear multifragmentation from heavy ion 
collisions involving asymmetric nuclei, it will be of interest to extend these dynamical models 
to include also the isospin degree of freedom and its fluctuations. 

This work was supported in part by the NSF Grant No. PHY-9509266. We would like 
to thank J.B. Natowitz for helpful discussions. 
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FIGURES 


Chemical & Mechanical Instability 



FIG. 1. Chemical (upper window) and mechanical (lower window) stability conditions as func¬ 
tions of density at fixed temperatures T = 4, 8 and 12 MeV for 5 = 0.0 (left panel), and 5 = 0.2 
(right panel). 
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Chemical & Mechanical Instability 



FIG. 2. Same as Fig. |] but for <5 = 0.4 (left panel), and 6 = 0.6 (right panel). 
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P (MeV/fm 3 ) 



FIG. 3. Pressure as a function of density at fixed temperatures T = 4, 8, 12 and 16 MeV for 
5 = 0.0, 0.2, 0.4 and 0.6. The isothermal (diffusive) spinodals are plotted using dashed (dotted) 
lines. 
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Chemical & Mechanical Instability 
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FIG. 4. Chemical (upper window) and mechanical (lower window) stability conditions as func¬ 
tions of 5 at a fixed temperature T = 10 MeV and various densities. 
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FIG. 5. Pressure as a function of density at a fixed temperature T = 10 MeV at various 6. The 
isothermal (diffusive) spinodals are plotted using dashed (dotted) lines. 
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p (fm 3 ) 6 


FIG. 6. Pressure as a function of density (left panel) and isospin asymmetry (right panel) 
at a fixed temperature T = 10 MeV along the boundary of diffusive spinodals (upper windows) 
and isothermal spinodals (lower windows) by using three forms of F(u ) and a compressibility of 
K = 380 MeV. 
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FIG. 7. Same as Fig. §j for a compressibility of K = 200 MeV. 
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